Towards a Bose-Glass of ultracold atoms in a disordered crystal of light 
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Starting from one- dimensional Mott Insulators, we use a bichromatic optical lattice to add con- 
trolled disorder to an ideal optical crystal where bosonic atoms are pinned by repulsive interactions. 
Increasing disorder, we observe a broadening of the Mott Insulator resonances and the transition to 
an insulating state with a flat density of excitations, suggesting the formation of a Bose-Glass. 
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Disorder plays a crucial role in statistical and con- 
densed matter physics and it contributes in a substantial 
way to the mechanism of transport and conduction. As 
originally predicted by Anderson disorder can lead 
to localization of a wave scattered by random impurities. 
Also interactions are well known to induce localization ef- 
fects, as happens in the Mott Insulator (MI) phase [21, |3|, 
in which a bosonic lattice system at zero temperature, in- 
stead of condensing in a superfluid (SF) state, forms an 
insulating "solid" with integer filling of the lattice sites. 
Much theoretical effort has been devoted to investigate 
the combined role of disorder and interactions in the SF- 
insulator transition observed in many condensed-matter 
systems [4]. If the disorder in the external potential is 
large enough, these systems are expected to enter an insu- 
lating state, the so-called Bose-Glass (BG), as predicted 
in the seminal paper Q , characterized by a gapless exci- 
tation spectrum and a finite compressibility [I, d, 0, S] • 

In this work we add controlled disorder to a collec- 
tion of ID ultracold ^''Rb gases in the MI phase by using 
a non-commensurate periodic potential superimposed on 
the main lattice, that introduces a randomization of the 
energy landscape on the same length scale as the lat- 
tice spacing. As a result, the characteristic resonances 
in the MI excitation spectrum are lost and the system 
rearranges to form a state with vanishing long-range co- 
herence and a broadband excitability. These observations 
constitute the first significative evidence in the direction 
of demonstrating a BG phase of ultracold atoms. 

At zero temperature the many-body quantum state of 
an interacting gas of identical bosons in a lattice poten- 
tial is well described by the Bose-Hubbard Hamiltonian 
H = -JJ2{j,r) % i^j - ^j^j^ where 

bj (b^j) is the destruction (creation) operator of one par- 
ticle in the j-th site, hj = V-hj is the number operator, 
and {j^j') indicates the sum on nearest neighbors [11,0]. 
The total energy results from the sum of three terms: J 
is the hopping energy^ proportional to the probability of 
tunnelling between adjacent sites, U is the interaction en- 
ergy^ arising from on-site interactions (repulsive for ^^Rb, 
for which /7 > 0), and Cj is a site-dependent energy ac- 
counting for inhomogeneous external potentials. 

When Cj = the ground state is determined by the 
competition between J and U. When U <^ J the system 
is in a SF state, in which the bosons are delocalized and 



tunnelling ensures long-range coherence. Instead, when 
^ J, the system is in a localized MI state, where phase 
coherence is lost and number Fock states are created at 
the lattice sites. The transition from a SF to a MI for ul- 
tracold atoms in an optical lattice has been first reported 
for a 3D system in Q and for an array of ID gases in [13]. 
The phase diagram of the system depends on the chemi- 
cal potential fi (related to the atomic density) and shows 
the existence of MI lobes with integer number of atoms 
per site (Fig. [1^). In the experiments, an additional 
harmonic confinement is present, resulting in a smooth 
variation of the density across the sample. As a result, 
in a single experimental run one typically averages on an 
extended range of local ja (vertical dotted line in Fig. [1^). 
This inhomogeneity precludes a sharp phase transition, 
due to the coexistence of SF and MI near the critical 
point, and is responsible for the formation of domains 
with different fillings deep in the MI phase [ill ]. 

In the presence of a disordered external potential an 
additional energy scale A enters the description of the 
system. We consider the case of bounded disorder, in 
which Cj G [—A/2, A/2]. In the presence of weak dis- 
order A < U the MI lobes should progressively shrink 
and a new BG phase should appear (Fig. [TJd), eventu- 
ally washing away the MI region for A>/7(Fig. [Ifc) 3. 
This BG phase shares some properties with the MI state, 
namely both are insulating states, with vanishing long- 
range coherence and vanishing superfluid fraction. How- 
ever, differently from the MI, the BG presents a gapless 
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FIG. 1: Phase diagram for disordered interacting bosons 0. 
Depending on the ratio between tunnelling energy J, interac- 
tion energy U and disorder A, the system forms a superfluid 
(SF), a Mott-Insulator (MI) or a Bose-Glass (BG) state. 
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FIG. 2: In a homogeneous MI the tunnelhng of one boson from 
a site to a neighboring one has an energy cost AE — U. In 
the presence of disorder this excitation energy becomes = 
U =b Aj, where Aj is the local site-to-site energy difference. 

excitation spectrum and a finite compressibility. 

In order to understand the physics happening when ap- 
proaching the BG phase, we consider the limit J ^ and 
unitary filling of the lattice sites. In a MI an energy gap 
exists, since the elementary excitation - the hopping of a 
boson from a site to a neighboring one - has an energy 
cost U. The presence of disorder introduces random en- 
ergy differences Aj G [—A, A] between neighboring sites 
(see Fig. ^jp) and the energy cost for such a process be- 
comes U ± Aj, that depends on the position. In the full 
BG, when A > [/, an infinite system can be excited at 
arbitrarily small energies and the energy gap shrinks to 
zero. Despite the zero energy gap, excitations only occur 
locally and the BG remains globally insulating. 

For a system of ultracold bosons in an optical lattice 
one can introduce controlled disorder in different ways. 
Experiments have already been realized with laser speck- 
les in combination with atomic Bose-Einstein conden- 
sates (BECs) pi,[i3,[li,[i5j. Optical disorder on much 
smaller length scales can be obtained by using two-color 
lattices 0, i-e. superimposing on the already existing 
lattice a second weaker lattice with non-commensurate 
spacing. 

In the experiment, an optical lattice {main lattice) is 
produced by using a Titanium: Sapphire laser at Ai = 830 
nm. Disorder is introduced by using an auxiliary lat- 
tice {disordering lattice) obtained from a fiber- amplified 
diode laser at A2 = 1076 nm. The resulting potential 
along the lattice axis x is V{x) = siEiiism^{kix) + 
52^i?2 sin^(/c2x), where si and 52 measure the height 
of the lattices in units of the recoil energies Ejn = 
h^/{2m\l) ^ hx 3.33 kHz and Er2 = h^/{2mXl) 
h X 1.98 kHz, h is the Planck constant and m the mass 
of a ^^Rb atom. When 52 <C si the disordering lattice 
has the only effect to scramble the energies e^-, which 
are non-periodically modulated at the length scale of the 
beating between the two lattices (2/Ai — 2/A2)~"^ = 1.8 
/im, corresponding to 4.3 lattice sites (see Fig. ^jp). In 
Fig. [St we plot a histogram of the site-to-site energy dif- 
ferences \Aj\ = \ej — ej-i \ occurring over the size of our 
samples (32 /im) for 52 = 1. Recent theoretical works 
have demonstrated that in finite-sized systems 
this quasiperiodic potential can produce the same effects 
induced by a truly random potential and allow the ob- 



servation of a BG of ultracold atoms. 

We perform the experiment with a collection of ~ 10^ 
independent ID bosonic systems tightly confined in the 
sites of a 2D optical lattice {trapping lattice) produced 
with laser light at Ai = 830 nm and aligned along y and 
z axes (see Fig. ^). The source of ultracold atoms is 
provided by a BEG of 1.5 x 10^ ^'^Rb atoms. In order 
to create the atom tubes we adiabatically increase the 
intensity of the 2D trapping lattice by using a 100 ms 
long exponential ramp with time constant 30 ms. After 
the ramp the lattice height, measured in units of Eri, 
is s± = 40 and the atoms are confined in the lattice 
sites with trapping frequency u± = 42 kHz along y and 
z. The confinement along the tubes is much more loose, 
being caused by the gaussian shape of the laser beams 
and the harmonic magnetic trap, resulting in a trapping 
frequency Ux = 75 Hz. The single-particle tunnelling rate 
between adjacent tubes is Jy = Jz = 0.4 Hz, that can be 
completely neglected on the timescale of the experiments. 

Together with the 2D lattice we switch on the bichro- 
matic lattice at wavelengths Ai and A2 along the direction 
of the tubes (see Fig. [3^) by using the same exponen- 
tial ramp. Then we characterize the many-body state 
of the ID systems by measuring the excitation spectrum 
and observing the interference pattern after time-of- flight 
(TOF), i.e. switching off the confining potentials and 
imaging the density distribution after expansion. 

The excitation spectrum of the ID gases in the bichro- 
matic lattice is measured with the Bragg spectroscopy 
technique introduced in [10]. A sinusoidal modulation of 
the main lattice height si with frequency u and ampli- 
tude 30% stimulates the resonant production of excita- 
tions with energy hv. We detect the excitations produced 
after 30 ms of modulation by first decreasing in 15 ms the 
intensity of the lattices back to si = = 5, S2 = in 
the 3D SF phase, waiting 5 ms and then switching off the 
potentials and imaging the atoms after TOF. The width 
of the central peak in the images is related to the energy 
transferred to the atomic system, giving information on 
the excitability of the sample at that frequency [lOj . 




FIG. 3: a) A collection of ID Rose gases in a deep 2D lattice 
experience the bichromatic potential V{x). b) Potential en- 
ergy along V{x) for si — 16 and S2 = 1. c) Histogram of the 
energy differences between neighboring sites for S2 = 1. 
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A typical spectrum for the MI state with si = 16 and 
52 = is shown in Fig. 2^. Here we detect the presence 
of an excitation peak at u = 1.9(1) kHz, corresponding to 
the interaction energy U. At twice the frequency 3.8(1) 
kHz we observe a second peak [3] , that can be attributed 
to higher-order processes and to excitations taking place 
at the boundary between different MI domains. For these 
parameters the peak chemical potential is ja 2.8/7 (see 
vertical line in Fig. [1]), giving rise to domains with 1, 2 
and 3 atoms per site. The same measurements are then 
repeated in the presence of the disordering lattice. Fig. 
[D3-e show the measured spectra for increasing disorder 
from S2 = 0.2 to S2 = 2.5. At S2 = 0.5 (Fig. Hb) one 
already detects the disappearance of the characteristic 
peak structure of the MI and the appearance of a broader 
spectrum. At the largest disorder height S2 = 2.5, when 
the maximum energy difference A between neighboring 
sites is 3.2 kHz c:^ 1.7U and one expects to have entered 
the full BG phase, the resonances are completely lost. 

Additional information can be obtained by analyzing 
the TOF images, that provide a measurement of phase 
coherence. We first prepare the system in a state with 
arbitrary si and 52, then we suddenly ramp in ^ 40 /is 
the lattice heights to si = 25 and S2 = (while keeping 
s± = 40) and finally switch off the confining potentials. 
This time is short enough not to change the coherence 
properties of the system, but allows us to project the 
state under investigation onto a same reference state [l^ . 
In Fig. [5^1 we report images of the density distribution 
after a TOF of 20 ms for 52 = 2.5 and different val- 
ues of si. The presence of vertical interference fringes 
is an indicator of long-range coherence along the tubes. 
When increasing the height of the main lattice, we ob- 
serve a progressive loss of coherence indicating the tran- 
sition from a SF to an insulating state [17]. In Fig. 
we report the coherent fraction, measured as the number 
of atoms in the interference peaks divided by the total 
number of atoms both for 52 = and for 52 = 2.5. 

The combination of the excitation spectra measure- 
ments and the TOF images indicates that, with increas- 
ing disorder, the system goes from a MI to a state with 
vanishing long range coherence and a flat density of exci- 
tations. The concurrence of these two properties cannot 
be found in either a SF or an ordered MI, and is consis- 
tent with the formation of a BG. Rigorously speaking, the 
BG phase should be characterized by a gapless spectrum. 
Detecting the absence of a gap is technically challenging, 
since it would require a measurement of excitability at 
arbitrarily small energies. A direct measurement of a 
small energy gap cannot be accomplished with the mod- 
ulation technique we have used, that works well only for 
frequencies u much larger than the reciprocal of the mod- 
ulation time r = 30 ms, i.e. starting from a few hundred 
Hz. However, the excitation spectrum is expected to be 
gapless only for an infinite system, while finite-sized sys- 
tems always have discrete energy spectra. Nevertheless, 
we expect the density of excitations in a finite-sized BG 
to lose the characteristic resonances of the MI and to be- 
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FIG. 4: a-e) Excitation spectra for si = 16 and different val- 
ues of S2. f) Effect of the modulation at frequency zy = 1.9 
kHz as a function of S2. The line in a) is a fit to the experi- 
mental points, while the curves in b-f) are calculated from a 
model of inhomogeneous broadening of the MI resonances. 



come flat [2] . In the following we will show that for weak 
disorder the spectra in Fig. [H can be explained with the 
inhomogeneous broadening of the MI peaks, accompa- 
nied by the consequent reduction of the gap, which is the 
first prerequisite for the formation of a BG. 

We have developed a model in which we calculate the 
inhomogeneous broadening of the MI resonances at U 
and 2U caused by the disordered distribution of site-to- 
site energy differences Aj (see Fig. [JJd). The broadened 
spectra are obtained by convolving the MI spectrum (fit- 
ted in Fig. [3^ with a double gaussian + a linear pedestal, 
solid line) with the distribution of energy shifts across the 
lattice. The result of this convolution is reported in the 
dashed curves of Fig. [4)3-e, showing a fine agreement 
with the experimental data for S2 ^ 1 [18.] . 

In order to more quantitatively analyze the disappear- 
ance of the MI resonances, we report in Fig. [3f the width 
of the density distribution after an excitation at u = 1.9 
kHz, corresponding to the first resonance in the MI spec- 
trum, as a function of 52 (empty circles). The solid line is 
the theoretical prediction based on the model of inhomo- 
geneous broadening, that nicely match the experimen- 
tal findings for 52 < 1.5. One expects this broadening 
to happen when approaching the transition from a MI 
to a BG, with the energy gap progressively closing with 
increasing disorder [6|. Eventually, when the broadened 
resonances reach zero-energy and the gap completely dis- 
appears, the transition to a BG is expected to occur. We 
find that the agreement with the model breaks down for 
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FIG. 5: a) Density distribution after time-of- flight for dif- 
ferent values of si and S2 = 2.5. b) Coherent fraction as a 
function of si for S2 = and S2 — 2.5 and sketch of the 
corresponding scans in the phase diagrams. Each point cor- 
responds to a different set of images 10 images per set). 



52 > 1.45 (vertical dotted line), when indeed the maxi- 
mum site-to-site energy difference is larger than U. At 
this point the atoms rearrange in a new insulating state 
with different sites filling, and the simple model of MI 
inhomogeneous broadening breaks down. 

Similar results for the excitation spectrum have been 
recently predicted in [16], where the authors study the 
response of a ID Bose gas to a periodic amplitude modu- 
lation of a superlattice. Our experimental findings are in 
good agreement with the numerical results of [l6| , show- 
ing a broadening of the MI resonances for weak disorder 
and their complete disappearance when entering the BG. 

Despite this remarkable change in the excitation spec- 
trum, an undeniable hallmark of the formation of a BG 
is still lacking. In particular, if a precise detection of the 
vanishing gap could be performed, it would discriminate 
between a genuine BG and a strongly inhomogeneous MI. 



To this aim a measurement of compressibility would pro- 
vide a decisive proof for the creation of such a state. 
However, in actual MI experiments the existence of mul- 
tiple domains leads itself to already globally compress- 
ible systems [11]. Regarding the coherence properties, 
few theoretical works predict that for restricted ranges 
of parameters one should observe an enhancement of co- 
herence due to the presence of disorder in the route to 
the BG pj]. As shown in Fig. [5)3, we do not observe any 
appreciable difference varying S2 within the reproducibil- 
ity limits of our experiment [20]. However, a quantitative 
prediction of these effects for real experimental systems 
is absent. Furthermore, in the progression from SF to 
MI/BG, one simultaneously probes extended regions of 
the phase diagram (as sketched in Fig. [SJd), that smooth 
the transition and could smear out these effects. 

In conclusion, we have added controlled disorder onto 
a collection of ID bosonic gases in an optical lattice by 
means of an additional lattice with different spacing. We 
have reported on the first observation of the transition 
from a Mott Insulator to a state with vanishing coher- 
ence and a flat density of excitations. These coexisting 
properties suggest the formation of a Bose-Glass, which 
is expected to appear for our parameters [tI, [sI, [l^] . Fu- 
ture work will be done in the direction of implementing 
new techniques for a more exhaustive characterization 
of this novel state, possibly including a measurement of 
compressibility and new detection schemes [21]. 
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